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f ABSTRACT 

The present paper is the effects of radiation, chemical reactions and thermo Diffusion (Soret 
effect) to MHD flow of an electrically conducting, incompressible, viscous fluid past an impulsively 
moving isothermal vertical plate through porous medium in the presence of uniform suction. Thermo 
diffusion is one of the mechanisms in transport phenomena in which molecules are transported in a 
multi component mixture driven by temperature gradients. A flow of this type represents a new class of 
boundary layer flow at a surface of finite length. The equations governing the flow field are solved by 
analytical method. The velocity, temperature, concentration and skin friction have been evaluated for 
variation in the different governing parameters. 
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I. INTRODUCTION 

Thermo diffusion along with molecular diffusion occurs in many engineering systems and in nature. 
Thermo diffusion has a great effect on the concentration distribution in a multi component mixture. The 
variations of composition and temperature may either lessen or enhance the separation in mixtures. The thermo 
diffusion phenomenon also plays a major role in the hydrodynamic instability analysis in mixtures, 
investigations of mineral migration, the mass transport modeling in living matters, and composition variation 
studies in hydrocarbon reservoirs. Diffusion is one of the major mechanisms of transport phenomena. Molecular 
diffusion is the movement of molecules from a higher to lower chemical potential. Thermo diffusion and 
pressure diffusion are additional ways in which molecules are transported in a multi component mixture driven 
by temperature and pressure gradients, respectively. The thermo diffusion phenomenon was discovered by 
Ludwig [1] and Soret [2], and named as the Soret effect. The Soret coefficient is the ratio of the thermo 
diffusion coefficient to the molecular diffusion coefficient. 

Furthermore, the study of thermo hyaline instability with thermo diffusion in a fluid saturated porous 
medium is of importance in geophysics, ground water hydrology, soil science, oil extraction. The reason is that 
the earth's crust is a porous medium saturated by a mixture of different types of fluids such as oil, water, gases 
and molten form of ores dissolved in fluids. Thermal gradients present between the interior and exterior of the 
earth's crust may help convection to set in. The thermal gradient in crude oil can have a strange effect on the 
distribution of petroleum components in an oil deposit. The thermal gradient causes the Soret effect which 
makes the larger molecule components to have a tendency to rise, while smaller molecule components go down 
to the bottom of the well. However, gravity causes the heavy components in a fluid to fall and the lighter ones to 
rise. This means that the distribution of components in a given well is neither consistent nor readily predictable. 
Pal and Talukdar [3] analyzed the combined effect of mixed convection with thermal radiation and chemical 
reaction on MHD flow of viscous and electrically conducting fluid past a vertical permeable surface embedded 
in a porous medium is analyzed. Mukhopadhyay [4] performed an analysis to investigate the effects of thermal 
radiation on unsteady mixed convection flow and heat transfer over a porous stretching surface in porous 
medium. Hayat et al. [5] analyzed a mathematical model in order to study the heat and mass transfer 
characteristics in mixed convection boundary layer flow about a linearly stretching vertical surface in a porous 
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medium filled with a viscoelastic fluid, by taking into account the diffusion thermo (Dufour) and thermal- 
diffusion (Soret) effects. Bhavana et.al [6]. The Soret effect on free convective unsteady MHD flow over a 
vertical plate with heat source. Chenna Kesavaiah et.al [7] Radiation and Thermo - Diffusion effects on mixed 
convective heat and mass transfer Flow of a viscous dissipated fluid over a vertical surface in the presence of 
chemical reaction with heat source. 

MHD flows has received the attention of several researchers due to its diverse applications in many 
branches of sciences and technology such as in planetary fluid dynamics, power generation systems, nuclear 
reactor thermal dynamics and electromagnetic materials processing. On the other hand, MHD flows with 
combined heat and mass transfer has numerous applications in chemical, mechanical and biological sciences. 
Some important application are stellar and solar structures cooling of nuclear reactors, interstellar matter, liquid 
metals fluid, power generation system, radio propagation and aerodynamics. 

The study of heat generation or absorption in moving fluid is important in problems dealing with 
dissociating fluids. Possible heat generation effect may alter the temperature distribution; consequently, the 
particle deposition rate in nuclear reactors, electronic chips and semi conductor wafers. Rahman et al. [8] 
studied thermo-micropolar fluid flow along a vertical permeable plate with uniform surface heat flux in the 
presence of heat generation. Ibrahim et al. [9] found the analytic solution of MHD mixed convection heat and 
mass transfer over an isothermal, inclined permeable stretching plate immersed in a uniform porous medium in 
the presence of chemical reaction, internal heating, Dufour effect and Hall effects. Pal and Chatterjee [10] 
studied heat and mass transfer in MHD non-Darcian flow of a micropolar fluid over a stretching sheet embedded 
in a porous media with non-uniform heat source and thermal radiation. Cortell [11] investigated theoretically the 
effect of viscous dissipation as well as radiation on the thermal boundary layer flows over non -linearly 
stretching sheets considering prescribed surface temperature and heat flux. Raptis and Perdikis [12] studied 
numerically the steady two dimensional flow of an incompressible viscous and electrically conducting fluid over 
a non-linear semi-infinite stretching in the Presence of a chemical reaction and under the influence of magnetic 
field. Rajagopal et al. [13] studied a Falkner-Skan flow field of a second-grade viscoelastic fluid. Shit and 
Haldar [14] investigated the effects of thermal radiation and Hall current on MHD free-con vective flow and 
mass transfer over a stretching sheet with variable viscosity in the Presence of heat generation/absorption. 

Therefore the objective of the present paper is to study the MHD heat transfer flow of a Newtonian 
fluid moving vertical surface with heat generation or absorption by considering the radiation effect along with 
Soret effect, uniform heat and mass flux and chemical reaction. 



We consider the steady, two-dimensional laminar, incompressible flow of a chemically reacting, 
viscous fluid on a continuously moving vertical surface in the presence of a uniform magnetic field, radiation 
absorption with heat generation, uniform heat and mass flux effects issuing a slot and moving with uniform 
velocity in a fluid at rest. Let the x-axis be taken along the direction of motion of the surface in the upward 
direction and y-axis is normal to the surface. The temperature and concentration levels near the surface are 
raised uniformly. The induced magnetic field, viscous dissipation is assumed to be neglected. Now, under the 
usual Boussinesq's approximation, the flow field is governed by the following equations. 



II. 



FORMULATION OF THE PROBLEM 
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Continuity equation 

du dv 

— +— =0 (1) 

dx dy 

Momentum equation 

2 2 

du du , d u aB n v 

u — +v — = g p (r'-rj)+ gp {c'-ci) + v — r - — u — — u (2) 

dx dy dy p K 

Energy equation 

f dT' dT ,s \ d 2 T' dq 
pC p \u +v \ = k ___!!_ Qo ( r »_ r j) (3) 

dx dy J dy dy 



8C dC d 2 C d 2 T' 



+ v — = D M — —-Kr'(C'-C:) + D T — (4) 



dx dy dy dy 

The relevant boundary conditions are 

u = u ,v - —v Q = const < 0 
dT q dC j" 

=--, = -— at y = 0 (5) 

dy k dy D 

u — > 0 , T — > , C — > C ^ as y — > co 

The radiative heat flux q r is given by equation (5) in the spirit of Cogly et.al [15] 

dq . 8 e. , 
-=4(T'-T')I where I=\K —dX, 

dy 0 " ST' 

K , - is the absorption coefficient at the wall and e - is Planck's function, / is absorption coefficient 

In order to write the governing equations and the boundary conditions the following non dimensional quantities 
are introduced. 

u yv n T'-T' C'-C ' K v 2 <tB\ vQ n 

u = ,Y = T = —,C = = - JL ^,M = °—,Q = — ^-y 

«„, v ( qv ~\ ( j" v } v 2 pv 2 pC v 2 



f qv ^ . ( j"v I 
v gp | | vgP I 

p MC " y Kr ' V n Q < J — r ^ kV ° J c V r V kv Q J rm 

P r = , a r = — ,Q i = ,Gr = , Sc = ,Gc = (o) 

k v o 9V 2 0 pC p u u x>l D M u u y 2 0 

D M K T j" Aqvl qD T 
u = , R = , S „ = 



Where u , v are the velocity along the x , y -axis, is constant obtained after integration conservation of mass in 
pre-non dimensional form not mentioned, v is the kinematic viscosity, g is the acceleration due to gravity, T ' 
is the temperature of the fluid, is the coefficient of volume expansion, C is the specific heat at constant 

pressure, is a constant, a is the Stefan-Boltzmann constant, T' exceeds the free steam temperature r ' , C is 
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the species concentration, is the wall temperature, C ' is the concentration at the plate, T ' is the free steam 
temperature far away from the plate, C[ r is the free steam concentration in fluid far away from the plate, Pr is 
the Prandtl number, G r is the Grashoff number, k is the thermal conductivity of the fluid, B 0 is uniform 
magnetic field strength, M is the magnetic field parameter which is the ratio of magnetic force to the inertial 
force. It is a measure of the effect of flow on the magnetic field. The termg 0 (T ' - rj ) is assumed to be the 

amount of heat generated or absorbed per unit volume. Q 0 is a constant, which may take on either positive or 
negative values. 

III. SOLUTION OF THE PROBLEM 

In view of (6) the equations (2), (3) and (4) are reduced to the following non-dimensional form 

d 2 U dU ( \ \ 

M+ — \U=-GrT-GcC (7) 

dY^ dY 



I M + — II 
I K ) 



d'T dT 

+ Pr -(Q + R)PtT=0 (8) 

dY 2 dY 

d 2 C dC d 2 T 

+ Sc KrScC + ScS =0 (9) 

2 o 2 

dY' dY dY' 

The corresponding boundary conditions can be written as 

dT dC 

U = 1, =-1, =-1 at Y = 0 

dY dY (10) 

U -> 0,T -> 0,C -> 0 as Y -> oo 

Where G r is the thermal Grashof number, Gc is the solutal Grashof number, Pr is the fluid Prandtl number, 
Sc is the Schmidt number K r is the chemical reaction, Q is Heat source parameter, M is the Magnetic 

parameter, S is the Soret number and R is the Radiation Parameter. 

The study of ordinary differential equations (7), (8) and (9) along with their initial and boundary conditions (10) 
have been solved by using the method of ordinary linear differential equations with constant coefficients. We 
get the following analytical solutions for the velocity, temperature and concentration. 

U = Jj e'" 1 ' + J 2 e'" 2 *' + J 3 e" u ~ + J 4 e'""" 

T = B l e"' ,y 

C = zy' 1 ' + Z^e'" 1 ' 

The computed solution for the velocity is valid at some distance from the slot, even though suction is applied 
from the slot onward. This is due to the assumption that velocity field is independent of the distance parallel to 
the plate. The fluids considered in this study are air (Pr = 0.71) and water (Pr = 7.0). 

Skin friction 

f dU ) 
dy 

= J , m 2 + J , m 2 + J 3 m 4 + / 4 m 6 
Nusselt number 
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Sherwood number 

(dC\ 

5/7=1 I = Z , m , + Z ,m 
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APPENDIX 
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IV. RESULTS AND DISCUSSION 

The variation of non-dimensional velocity U for different values of Grashof number (Gr), Solutal 

Grashof number ( G c ) , Porosity permeability ( K ) , Chemical reaction parameter ( K r ) , Magnetic parameter 
( M ) , Prandtl number ( P r ) , Heat source parameter ( Q ) , Radiation parameter ( R ) Soret number ( 5 0 ) and 
Schmidt number (Sc) are shown in figures (1) - (10). 



Figure (1) shows the influence of thermal buoyancy force parameter ( G r ) on the velocity. As can be 

seen from this figure, the velocity profile increases with increases in the values of the thermal buoyancy. We 
actually observe that the velocity overshoot in the boundary layer region. Buoyancy force acts like a favourable 
pressure gradient which accelerates the fluid within the boundary layer therefore the solutal buoyancy force 

parameter ( G c) has the same effect on the velocity shown in figure (2) asG r . Figure (3) shows that the axial 

velocity increases with a rise in permeability parameter (AT ) . Moreover, for small permeability, the axial 

velocity reduces faster than when the permeability is higher. Figure (4) depict effect of chemical reaction on the 
flow variables is shown. In figure (4) velocity is clearly boosted with stronger chemical reaction i.e. as the 

chemical reaction parameter ( K r ) increases profiles are lifted continuously throughout the boundary layer, 

transverse to the plate. A distinct velocity escalation occurs near the wall after which profiles decay smoothly to 
the stationary value in the free stream. Chemical reaction therefore boosts momentum transfer i.e. accelerates 

the flow. Figure (5) reveals that axial velocity decreases as Magnetic parameter ( M ) increases. This may be 

attributed to the fact that an increase in M signifies an enhancement of Lorentz force, thereby reducing the 
magnitude of the velocity. This figure further indicates that blood velocity in the capillary decreases with 

increase in distance from the lower wall of the capillary. Figure (6) illusttates the effect of Prandtl number ( P r ) 

on the velocity. It is noticed that as the Prandtl number ( P r ) increases, the velocity decreases. As seen in the 
earlier cases, far away from the plate, the effect is much significant. Figure (7) illustrate the variation of velocity 
with the effects of heat absorption parameter ( Q ) . The velocity values are clearly reduced with increasing Q ; 
again an overshoot is computed close to the plate both in the presence and absence of heat source. Heat source 
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however suppresses the overshoot. Figure (8) show the effect of radiation parameter ( R ) on the velocity it is 

observed that the velocity decrease as the radiation parameter ( R ) increases. This result qualitatively agrees 

with expectations. Figure (9) shows the variation of the velocity boundary-layer with the Soret number ( S r ) .It 

is found that the velocity boundary layer thickness increases with an increase in the Soret number. Figure (10) 

illustrates the effect of the Schmidt number (Sc) on the velocity. The Schmidt number (Sc) embodies the ratio 

of the momentum diffusivity to the mass (species) diffusivity. It physically relates the relative thickness of the 
hydrodynamic boundary layer and mass-transfer (concentration) boundary layer. It is noticed that as Schmidt 

number ( Sc) increases the velocity field decreases. 

Figures (11) - (13) give some characteristic temperature profiles for different values of Prandtl 
number ( Pr ) , heat source parameter ( Q ) and radiation parameter ( R ) . Figure (11) presents the change in the 

temperature distribution in the boundary layer, when the Prandtl number ( P r ) changes gradually. It shows that 

as the Prandtl number increases, the temperature of the boundary layer diminishes. This may be attributed to the 
fact that the thermal boundary layer thickness reduces with an increase in Prandtl number. Further, this figure 
indicates that the temperature gradient at the surface increases with a rise in Prandtl number. This implies that an 
increase in Prandtl number is accompanied by an enhancement of the heat transfer rate at the wall of the blood 
vessel. The underlying physics behind this can be described as follows. When blood attains a higher Prandtl 
number, its thermal conductivity is lowered down and so its heat conduction capacity diminishes. Thereby the 
thermal boundary layer thickness gets reduced. As a consequence, the heat transfer rate at the vessel wall is 

increased. Figure (12) illustrates the effect of the heat source parameter ( Q ) on the temperature. It is noticed 

that as the heat source parameter increases, the temperature decreases. Figure (13) illuminates a very important 
effect of thermal radiation on the temperature profile. This figure emphasizes that as thermal radiation increases 
during blood flow in capillaries; there is a significant rise in the thickness of boundary layer. Thereby the 
temperature of the boundary layer in enhanced by an appreciable extent. 

The effects of order of chemical reaction on concentration distributions are displayed in figure (14) 
respectively. It is observed from these figures that an increase in the chemical reaction order n leads to decrease 
the concentration profiles. Figure (15) illustrate the effect of the Prandtl number on the concentration 
profiles respectively. It is evident that increases in the Prandtl number leads to increases on the 
concentration profiles. 

From figure (16) Soret number (S 0 ) increases with increasing values on the concentration profiles. 
Figure (17) illustrate the concentration field distributions with transverse coordinate for different Schmidt 
number( 5c ) . An increase in Sc causes a considerable figure (17). A much greater reduction is observed in 

concentration values in figure (17). An increase in (Sc) will suppress concentration in the boundary layer 
regime. 

Higher Sc will imply a decrease of molecular diffusivity D causing a reduction in concentration 
boundary layer thickness. Lower Sc will result in higher concentrations i.e. greater molecular (species) 
diffusivity causing an increase in concentration boundary layer thickness. For the highest value of Sc , the 
concentration boundary layer thicknesses are of the same value approximately i.e. both species and momentum 
will diffuse at the same rate in the boundary layer. 

The effects of Soret number ( S 0 ) versus Grashof number (G r) on the local skin-friction coefficient shown 
in figure (18). From these figure we observe that the local skin-friction coefficient increase as Soret number 
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CONCLUSIONS 

The problem of steady laminar, hydromagnetic two-dimensional mixed convection flow due to stretching 
sheet in the presence of Soret effects was investigated. The resulting partial differential equations, describing the 
problem, are transformed into ordinary differential equations by using similarity transformations. These 
equations are more conveniently solved analytically by perturbation technique for the computation of flow, heat 
and mass transfer characteristics, skin-friction coefficient, for various values of the magnetic field parameter, 
heat source, Soret numbers, etc. Comparisons of the present results with previously published work on some 
limiting cases were performed and results were found to be in excellent agreement. From the present 
investigation following conclusions are drawn: 

• Increasing the magnetic field parameter decreases the flow velocity near the stretching sheet in the 
momentum boundary layer, whereas its effect is reversed as the space variable of porous permeability 
parameter. 

• Temperature deceases where as increases in radiation parameter. 

• Skin friction increases with increase in the value of Soret number. 

• Temperature decreases with increase in non-uniform heat source parameter 

• Increase in the value of Schmidt number the concentration profiles decreases. 

It is hoped that the present investigation is useful in the study of MHD flow over a vertical plate embedded in 
non-Darcian porous medium which can be utilized as the basis for many scientific and industrial applications 
and studying more complex problems involving Soret effects. 
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Figure (1): Velocity profiles for different values of Gr 
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Figure (2): Velocity profiles for different values of Gc 
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Figure (3): Velocity profiles for different values of K 
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Figure (4): Velocity profiles for different values of Kr 
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Figure (5): Velocity profiles for different values of M 
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Figure (6): Velocity profiles for different values of Pr 
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Figure (7): Velocity profiles for different values of Q 




Figure (8): Velocity profiles for different values of R 
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Figure (9): Velocity profiles for different values of S g 
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Figure (11): Temperature profiles for different values of Pr 
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Figure (12): Temperature profiles for different values of Q 
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